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EXAMPLE FOR COMBUSTION PRODUCTS OF BORON-CONTAINING FUEL 

By Leonard K. Tower 


SUMMARY 

Jet engine fuels containing boron form boric oxide (B 2 O 3 ) when 
burned. This substance can exist both as a liquid and a vapor at high 
temperatures. If a portion of the BgOj condenses during expansion, the 
performance of the engine can be affected appreciably. The physically 
possible effect of condensation on performance is often bracketed by 
assuming that expansion occurs either with no change in composition or 
with equilibrium condensation. 

One method of calculating the isentropic expansion process is by 
means of equations relating pressure, temperature, and density. Commonly 
used equations from engineering thermodynamics relating these variables 
are not suited to the analysis of an isentropic process involving equi- 
librium condensation. This report presents equations relating tempera- 
ture, pressure, and. density for an isentropic process with equilibrium 
condensation. 

The importance of using equations which properly describe this proc- 
ess is emphasized by a calculation discussed herein of the isentropic 
variation of pressure with temperature involving combustion products of 
boron. At the chosen condition, the coefficient for the isentropic vari- 
ation of pressure with temperature was 5.8 with no phase change assumed 
and 12.9 with condensation assumed to occur. 

A procedure is given for integrating the equations involving pres- 
sure, temperature, and density in order to analyze an expansion process. 
An example in the form of a problem concerning the use of ethyldecaborane 
in a ramjet engine serves to illustrate this procedure. 


INTRODUCTION 

The large heating values of boron-containing fuels offer the possi- 
bility of improvement in the flight range of air-breathing jet engines. 
The degree of improvement over hydrocarbon fuels depends in part both on 
the thermal properties of the combustion products and on the nature of 
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the expansion process. The "boron-containing fuels differ from conven- 
tional fuels "because they produce "boric oxide (B 2 O 3 ), which can exist 
simultaneously as a vapor and a liquid. If condensation of the oxide 
occurs during the expansion process, an appreciable effect upon thrust 
may be encountered under certain conditions. 

The degree to which condensation occurs in actual expansion processes 
depends upon the rate of condensation relative to the time of expansion. 

In the absence of information on condensation rates it is customary to 
bracket the region of probable performance with calculations that assume 
either infinitely fast or infinitely slow condensation rates. These are 
designated respectively as expansion with equilibrium condensation and 
expansion with no condensation. 

Expansions with no condensation are readily handled and understood 
by means of elementary thermodynamic concepts. Expansions with equilib- 
rium condensation are more difficult to handle. It is important that 
performance calculations involving the assumption of phase equilibrium 
be based upon correct equations and procedures . 

The purpose of this report is to define certain properties for an 
expansion process involving combustion products with equilibrium conden- 
sation. These properties are (l) the variation of pressure with density 
at constant entropy (dp/dp) ( 2 ) the variation of pressure with tempera- 
ture at constant entropy (dP/dT) s j (3) the variation of stream velocity 
at constant entropy (dU) s . A numerical integration of these terms pro- 
vides the pressure, temperature, density, velocity, and local Mach number 
at any point during the expansion. A procedure for integrating these 
terms is given, together with an example of the analysis of a specific 
expansion condition. 

The equations and procedure herein omit the effect of dissociation 
of gases. Dissociative reactions generally assume importance only above 
4000° R at pressures encountered in high-speed flight. For practically 
all conditions where boron-carbon-hydrogen fuels might be burned in air- 
breathing jet engines, saturation of the gases in the exhaust nozzle will 
not occur until a temperature less than 4000° R occurs. 

Previous reports present other methods whereby the expansion process 
with phase equilibrium can be analyzed. Reference 1 contains a general 
procedure by which the nozzle outlet velocity can be obtained for speci- 
fied nozzle inlet conditions and expansion ratio. Isentropic expansion 
can be considered for systems with compositions either frozen or in phase 
and chemical equilibrium. In addition, reference 1 presents a procedure 
for determining combustion temperatures which includes the effect of 
dissociation. 
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Because of the effort required in applying the procedures of refer- 
ence 1 to specific performance calculations, the authors of reference 2 
developed a simplified method for fuels containing boron, carbon, and 
hydrogen burned in air. The thermal properties for combustion products 
of certain stoichiometrically burned fuels are added to obtain properties 
for fuels containing boron, carbon, and hydrogen in any proportion and at 
any equivalence ratio. The properties are in tabular form and also in 
the form of pressure-enthalpy charts. 

The charts of reference 2 include the effect of both vaporization 
and phase change. In regions where two phases are present, the charts 
can properly be used for the specific fuels and equivalence ratios for 
which they were calculated. The tables of reference 2 neglect dissocia- 
tion of gases, but they can be used to include the effect of phase change. 
It is correct to add properties tabulated for combustion products of 
specific fuels to obtain the properties for any boron- carbon-hydrogen 
fuel when equilibrium phase change occurs . However, the determination 
of conditions at many points in the expansion process involves numerous 
trial and error calculations. Furthermore, the Mach number cannot be 
determined in a manner consistent with the assumed process of equilibrium 
condensation from data contained in the tables. 

The equations and procedure in the present report provide a method 
alternative to those of references 1 and 2 for calculating the expansion 
process with equilibrium condensation of B2O3. In analyses where many 
points during the expansion are desired, this method may prove more con- 
venient than the others. Moreover, the equations presented herein greatly 
aid in visualizing the effect of condensation on the expansion process. 


DISCUSSION 


Since the purpose of this report is to define thermal properties 
useful in analyzing an isentropic expansion process with two phases in 
equilibrium, these properties will first be listed. The isentropic vari- 
ation of pressure with density is defined as 



( 1 ) 


where 



(2) 
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and 


e 


fr = 



( 3 ) 


Symbols are defined in appendix A. The definitions of Xy and the last 
expression are described more explicitly in appendix B. The isentropic 
variation of pressure with temperature is 



( 4 ) 


where 


e = 


'fr 


1 + 


1 - x v \rt ) 
Xy /AH y \~ 
1 - Xy \RT ) 


( 5 ) 


The differential form of the general energy equation for isentropic 
expansion can be expressed as 


( UdU ) s = - if dT = " if 1 d ( ln p ) (e) 

The speed of sound is 

c = VrP/p (7) 

and Mach number is merely 

M = U/c (8) 

Derivations of equations (2) , (5) , and (6) are collected in appendix C. 

Equations (2) and (5) for y and e y respectively, are rigorous 
where condensation is the only change in composition. If the composition 
is assumed to be frozen during expansion so that no phase change occurs, 
those terms in equations (2) and (5) involving the heat of vaporization 
AHy are zero. Equations (2) and (5) become the familiar expressions 


r r fr - ( c p/ c y) fr 


( 9 ) 
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e 



from which 


( 10 ) 



(ID 


is derived. 

When values of y and e computed for equilibrium phase change are 
compared to Xf r and for frozen composition, the marked effect of 

condensation on the flow process is emphasized. As an extreme example, 
consider the combustion products of boron at an equivalence ratio of 0.9, 
a pressure of 1 atmosphere, and a temperature of 3960° R (not the flame 
temperature) . Assume that the combustion products are not dissociated 
at this temperature but that phase equilibrium exists. The ratio of 
gaseous to liquid B 2 0 3 is then 2.102 at this condition. If expansion is 

assumed to occur with no condensation (frozen composition), Xf r and 

from equations (9) and (10) are 1.2083 and 5.802, respectively. If equi- 
librium condensation is assumed during expansion, x and 6 from equa- 
tions (2) and (5) are 1.1234 and 12.939, respectively. 

Values of y and e calculated from equations (2) and (5) include 
the effect of condensation, while the presence of dissociated gases is 
ignored. By means of the method presented in reference 1, both phase 
change and recombination of dissociated gases can be evaluated, although 
greater effort is involved than in the use of equations (2) and (5). For 
the example under consideration, y is 1.1214 and e is 13.091 if both 
phase change and recombination of dissociated products are included. 

These values agree reasonably well with the values neglecting dissociation 
calculated by use of equations (2) and (5). For the combustion products 
of most boron-containing fuels, saturated mixtures of B 2 0 3 vapor to which 
equations (2) and (5) are applicable can exist only under conditions where 
dissociation is small. 

By means of equations (2) , (5) , and (6) , expansion processes involv- 
ing equilibrium phase change can now be analyzed with the assumption that 
gaseous products are not dissociated. For small intervals of the expan- 
sion process, pressure and temperature do not change greatly. Both y 
and s can then be treated as constants in the interval. Equations (l), 
(4), and (6) can then be integrated for this small interval to give these 
equations : 


p z /p y = (p 2 / Py y 


(12) 


6 


MCA RM E57C11 


P z/ P y = (T z /T y ) 6 


(13) 



(14) 


These computations are repeated for the number of intervals into which 
the expansion is divided, with y and e recalculated for each interval. 


An example is included to illustrate the manner in which the equa- 
tions developed herein are used to analyze a one -dimensional expansion 
with phase equilibrium. The example shows the effect which the assump- 
tion of equilibrium phase change has on the expansion process. 

Consider the following problem: a ramjet engine burning ethyldeca- 

borane fuel operates at 60,000 feet, a flight Mach number of 4.0, and an 
equivalence ratio of 0.6. Inlet and diffuser pressure recovery is 0.443, 
and combustion efficiency is 1.00. Combustor inlet Mach number is 0.175. 
Flameholder drag is neglected. The exhaust-nozzle outlet temperature and 
velocity are desired for isentropic expansion to ambient pressure with 
equilibrium condensation. Also desired are nozzle contours required for 

(1) constant change of velocity with distance along the nozzle and 

(2) constant change of Mach number with distance along the nozzle. 

Charts and procedure discussed in appendix D were used to determine 
conditions at the combustor outlet. The combustion temperature and pres- 
sure were found to be 4400° R and 4.524 atmospheres, respectively 
(l atm = 2116 lb/sq ft). Dissociation of gases was neglected in this 
calculation. Since no condensed B 2 O 3 was present at the combustor outlet, 

the condition in the nozzle at which condensation became possible was 
found by a procedure in appendix D. For convenience in discussion, this 
condition where condensation can theoretically begin is referred to as 
nozzle station 1. Subsequent stations for the stepwise calculation in 
the condensed region are referred to as 2, 3, 4, and so forth. Eight 
steps were used to calculate that portion of the expansion process after 
the theoretical saturation point. The nozzle inlet was designated as 
station C, the nozzle outlet as station D. Nozzle station 9 thus corre- 
sponds to the nozzle outlet, station D. 

Values for P, T, U, c, M, y, and e for station 1 onward are pre- 
sented in table I. The detailed calculations shown in table II are dis- 
cussed more fully in appendix D. By following the headings of table II 
the reader can make similar calculations for other conditions and other 
fuels containing boron, carbon, and hydrogen. 


EXAMPLE 
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Observe in "table I that two stations, 1 and l(a), are shown to have 
the same condition of static temperature and pressure. The temperature 
at station 1 was assumed to be infinitesimally above the temperature 
where condensation starts in this particular case. Tha t at station l(a) 
was assume to be infinitesimally below the condensation condition. Sharp 
discontinuities in y and £ result because of the sudden appearance of 
terms containing ABy in equations (2) and (5) with the onset of conden- 
sation. The stream properties P, T, and U are the same at stations 1 
and 1(a) . However, the sonic velocity and Mach number jump because of 
the discontinuity in y with the onset of condensation. 

Between stations 1 and l(a) there is a 100-percent increase in £ 
for equilibrium condensation over the value for no condensation at the 
same temperatures. There is an 8-percent decrease in y for the same 
circumstances . 

From the data of table I were computed the nozzle contours required 
for (l) constant change of velocity with distance along the nozzle axis 
and (2) constant change of Mach number with distance. The nozzle was 
assumed to be of circular cross section and 2 feet long. The nozzle 
radius required to pass 1 pound of mixture per second is shown plotted 
against axial distance for both cases in figure 1. Also shown in figure 
1 is the nozzle contour for constant increase of velocity with distance 
if no condensation were to occur. 

For the case of constant Mach number increase a procedure in appendix 
D avoids a discontinuity in the nozzle contour between stations 1 and 1(a). 
Such a discontinuity would otherwise result from the sudden increase in 
r caused by the onset of condensation. Between the nozzle inlet (station 
C) and station 1 the value of e^> r used in computing stream velocity 
from equation (14) is assumed to be 5.15, while Tffj. is assumed to be 
1-24 in equation (7) for the speed of sound. The flow areas and Mach 
numbers between station C and the saturation condition calculated using 
these values of y^ r and e^ r assume no dissociation of gases. Incor- 
porating the effect of dissociation between these stations changes the 
flow areas and Mach numbers very slightly in this particular example. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, March 15, 1957 
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APPENDIX A 
SYMBOLS 

A specific nozzle area, sq ft/(lb/sec) 
a constant in expression for molar entropy 

Cp constant-pressure specific heat of substance, Btu/(lb mole) (°R) 
c local speed of sound, ft/sec 

Cp constant -pres sure specific heat of combustion products at constant 

composition, Btu/(lb)(°R) 

cy constant -volume specific heat of combustion products at constant 
composition, Btu/ (lb) (°R) 

g gravitational constant, 32.17 ft/sec^ 

Hp total enthalpy of combustion-product constituent (sum of chemical 
energy and sensible enthalpy), Btu/lb mole 

AHy heat of vaporization, Btu/lb mole 

h total enthalpy (sum of chemical energy and sensible enthalpy), 
Btu/lb 

Ahj, heat of combustion of fuel, Btu/lb 

J conversion factor, 778.16 ft-lb/Btu 

7 distance along nozzle axis, ft 

M Mach number, ratio of stream velocity to local speed of sound 
m molecular weight, lb/lb mole 

1 % mean molecular weight of combustion products, lb/lb mole 
n number of moles 

n’ number of moles of combustion-product constituent resulting from 
burning fuel stoichiometrically in 1 pound of air 

P static pressure of combustion products, atm (l atm = 2116 lb/sq ft) 
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p partial pressure of combustion-product constituent, atm 
R gas constant, 1.98718 Btu/(lb mole) (°R) 

S entropy of substance, Btu/(lb mole) (°R) 

s entropy of mixture of materials, Btu/°R 

T static temperature, °R 

U axial stream velocity, ft/sec 

V volume of mixture of combustion products, cu ft 
W weight, lb 

X number of moles of constituent divided by total moles of gas and 
vapor 

X coefficient for variation of pressure with density at constant 
entropy 

e coefficient for variation of pressure with temperature at constant 
entropy 

t arithmetic average of e over temperature interval 

p density, lb/cu ft 

$ equivalence ratio, actual fuel-air ratio divided by stoichiometric 

fuel-air ratio 

Subscripts: 

a air 

B combustor inlet 

C combustor outlet or nozzle inlet 

c condensed phase 

D nozzle outlet 

f fuel 

fr no change in composition (frozen composition) 
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g 

1 

k 

P 

s 

st 

T 

v 

y,z 

a,3,X 
T)=l, 2,3 


noncondensing gases only, vapor of condensing material excluded 

all gases, including vapor of condensing material 

all materials, including condensed phase 

constant pressure 

constant entropy 

stoichiometric 

. - ' < 

constant temperature 

vapor phase of condensing material 

adjacent stations in flow 

atoms of boron, hydrogen, and carbon, respectively, in fuel 
formula B a HpC\ 

stations in nozzle after condensation begins 
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APPENDIX B 


RELATIONS FOR ONE-DIMENSIONAL ISENTROPIC EXPANSION WITH 
NO CONDENSATION AND NO DISSOCIATION 

This appendix discusses the definitions of y and e for the isen- 
tropic expansion of two-phase systems with no condensation. Similar 
material has "been presented elsewhere (ref. 3), hut the topic is dis- 
cussed here because some of the terms derived are required in the analy- 
sis of equilibrium condensation in appendix C. 

Consider a mixture of ideal gases containing finely divided condensed 
material too large for Brownian motion. The ideal gas law applied to the 
gases in the mixture is 


PV = RJT2n i 


(Bl) 


The density of the mixture including the condensed phase is 

p = Sn^/V (B2) 

provided the volume of the condensed phase is negligible compared to that 
of the gas. Substituting equation (B2) in equation (Bl) gives 


P 2ni 

- = R JT * — 

p 2n k m k 


(B3) 


The term Z^mj^/Zn^ can be regarded as n^, the mean molecular weight of 

the mixture of gases and condensed material, as contrasted with the mean 
molecular weight of the gases alone, which is defined as Zr^m^/Zn^. 

The derivative (5 p/St) s is found by considering the general expres- 
sion for reversible entropy change 


T ds = dh - \ dP (B4) 

J 

If composition is fixed (no condensation) and there is thermal equilibrium 
between particles and gas, 

dh = Zn k (C°) k dT (B5) 

Set ds equal to zero. Then combine equations (B3), (B4), and (B5) to 
obtain 
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s,fr 


P 2n k( g p)k 
T R2n, 


(B6) 


Equation (B6) can "be expressed somewhat differently. Let the mole 
fraction for all constituents including the condensed phase he defined as 
X k = n k /2n i . The number of moles of condensed material is not included in 

the summation Znp. For this reason, ZX k = 1 + X c . This convention, 
used throughout the remainder of this report, gives 


zx k (cg) k 


s,fr 


l _ e l 

T ~ f r t 


(B7) 


The derivative (Sp/Sp) s can be found by differentiating the loga- 
rithm of equation (B3) with respect to density at constant entropy with 
2 n k m k constant for a given weight of mixture: 


/Sp\ p p /St\ /dp\ ^ p /^ 2n i\ 

Ws = p * Ws \H 2% V^T/ S 


(B8) 


For frozen composition (dZnp/dpJg is zero. Combining equations (B6) and 
(B8) gives 


/sr\ p 

2n k (cg) k 

(^)s,fr = rf P ‘ 

2n k(cg)k - *2% 


(B9) 


Thus, for frozen composition 


2n k ( C p) k/ 2n k m k 


** ~ t2n k (C-) k /Zn k m k ] - RZ^/Zn^ 


(BIO) 


The numerator of expression (BIO) is Cp and the denominator is cy for 
a unit weight of the mixture of gases and condensed material. The defini- 
tion of x for the case of frozen composition is thus extended by equa- 
tion (BIO) to include the situation of small particles in velocity and 
temperature equilibrium with the gas. 


Equation (BIO) expressed in terms of mole fraction is 

zx k (Cp) k 


Tfr = 




R 


(BH) 
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Equations (B7) and (Bll) can be combined to give 


fr 


Tfr 

T fr - 1 


(B12) 


Only by specifying frozen composition is the relation between s and y 
obtained as given by equation (B12) . 
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APPENDIX C 

RELATIONS FOR ONE-DIMENSIOML ISENTROPIC EXPANSION WITH 

PHASE CHANGE OF ONE CONSTITUENT AND NO DISSOCIATION 

In this appendix, equations (2), (5), and (6), for r, e, and dU, 
respectively, are derived. These expressions pertain to two-phase isen- 
tropic expansion with phase equilibrium. 


Derivation of e 


It is convenient to derive e first. The result can then be used 
to find x • Eo* 1 a fixed mass of material with gases behaving ideally, 

s = (Cl) 

where the summation over k terms includes the condensed phase and its 
vapor as well as the noncondensing gases. 


For the noncondensing gases and the vapor phase of the condensing 
material, collectively represented by the subscript i. 



dT - R In p^ + a^ 


(C2) 


For the condensed phase the entropy can be considered nearly independent 
of pressure. Hence (see ref. 4), 


f( C p)c 

S c =J^~ dT + a c 


(C3) 


In an isentropic process 


ds = 0 = Zn^dS^ + ZS^dn^ 


(C4) 


Representing both and n^ as functions of the temperature and 

pressure of the system yields 
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and 



Combining equations (C4) , (C5) , and (C6) gives 


(C6) 



(C7 ) 


If equation (C7) is to be valid, the particles must be in thermal equi- 
librium with the gas and vapor. 


The derivatives are now evaluated. First find Zn k (dS k /dp) T . From 
equations (C2) and (C3) 



(C8) 


and (SS c /Sp)j = 0. The partial pressure of a constituent in a mixture 
of ideal gases is 

Pi = Pn./Zni (C9) 


where 


Z^ = + Zn g (CIO) 

Assume that neither ng nor Zng changes for the noncondensing gases. 
Differentiate the logarithm of equation (C9) for the vapor phase and 
combine the result with equation (CIO) to obtain 



However, the partial pressure of a saturated vapor behaving as an ideal 
gas is a function of temperature only. The relation is known as the 
Clausius -Clapeyron equation (see ref. 4). In reference 4 it is shown 
that (dp v /dP) T -+0 if the molar volume of the liquid becomes negligible 
compared to that of its vapor. Equation (Cll) then gives 
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dn v\ _ 

5p“/ t “ P(l - Xy) 


(C12) 


where X v = n v/^ n l • For "the no nconde using gases, differentiation of 
equation (C9) and substitution of equation (C12) yield 


1 /^Pg\ 1 x v 1 

Pg W7 t " p + P U - V " p il - Xy) 


(C13) 


Finally, since (dp v /dP) T is zero, 



R Zn K 


Substituting for Zn from equation (CIO) results in 

o 


2n k 




(C14) 


The term 2 (c)n k /3p) ip of equation (C7) is considered next. The 
mole numbers of only the condensed phase and its vapor are assumed to 
change. The relation between moles of liquid and moles of vapor is 
dny = - dn c . Then 


2s k 



- - s => (sf) t 


(015) 


But 


Sy - s c = ABy/T (C16) 

Equation (C16), together with equation (C12), when substituted in equa- 
tion (C17) gives 



n v^v 
PT(l - Xy) 


which can be written 


A 

2S kb 

\\ Rny(AHy/RT) 

? / T P (1 • Xy) 


(C17) 
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The term ZS^dn^/d^p of equation (C7) is now considered. Differenti- 
ate the logarithm of equation (C9) for the vapor phase with respect to 
temperature : 

_l/ S Pv\ _l/^ n v\ _ 1 /dn v \ 

VyW ' /p W /p " 2 n iW ~. /P 

This yields finally 

Sn v\ n v 

St“/ p " p v (i - x v )V5t"/ p 

For a liquid surface in equilibrium with its vapor, the Claus ius-Clapeyron 
equation is 


(CIS) 


(C19) 


^Pv 

dT~ 


PvAHy 

RT 2 


(C20) 


if the molar volume of the liquid is negligible compared to that of its 
vapor (see ref. 4). This expression is applicable to vapors behaving as 
ideal gases. If the radii of the liquid surfaces are not too small, the 
value of AH V can be taken as that of a flat liquid surface. For in- 
stance, the surface energy of a drop as small as 2.2x10“® centimeter is 
negligible (ref. 5) . The ratio of vapor pressures of the small drop and 
the flat liquid is a function of surface energy of the small drop (ref. 
4). Therefore particles as small as 2.2x10"® centimeter have nearly the 
same vapor pressure as the flat surface. Substituting equation (C20) in 
equation (C19) results in 

/Sn v \ _ n v (AH v /RT) 

\3t~ / p ■ T(1 - S3 

Because only n v and n c change, 

LSk (^T L ) p = ^ 


( 021 ) 


(C22) 


Substitution of equations (C16) and (C2l) in equation (C22) gives 


HoT / p T(l - x^) 


(C23) 


The term Znj c (SSj J /^T)p of equation (C7) can now be found. From 
equations (C2) and (C3) , 
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and 

(ir)p = (c p } c/ t 

for any constituent. From equation (C9) 



(C24) 


for noncondensing gases. This expression, together with equation (C18) 
and equations (C24 ) , upon summation yields 


2n k(dT )p Snjj^pJ^/T (C25) 

Substituting equations (C14), (C17), (C23), and (C25) in equation (C7) 
yields 


2 n k( C p)k ^(AH^/RT ) 2 
/dp\ = T + T(1 - Xy) 

\ St /c "r v _" . Bn v (AH v /RT) 

P 2 n i + P(1 - Xy.) 


Let X k - n^/Zn.^ 
particles; then 


for any constituent including the condensed 



P 

T 


1 + (AH^/ET) 


(C26) 


In appendix B it is shown that if the composition of a mixture of gases 
and condensed particles is frozen (no change in composition) 


s 


fr 



ZX k (C°) k /R 
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Rewrite equation (C26) as 


e 



T 

P 


'fr 


1 + 


Xy 

1 - X v 
Xy 

1 - Xy 


(AH V /RT) 2 
(AHy/RT) 


( 5 ) 


Observe that e fr in equation (5) must be evaluated by insertion of those 
values of the X k ! s required for phase equilibrium at the temperature 
and pressure in question. 

Equation (5) was derived on the assumption that the mixture was 
saturated so that BgOj condensed. If the mixture is superheated or if 
expansion is frozen, no B 2 O 3 condenses, and the heat of vaporization is 
not recovered. Equation (10) in the DISCUSSION results. 


Derivation of y 

Consider equation (B 8 ), the differential form of the ideal gas law: 

(!) s ”? + i (l) s + ^(^r) s (B8) 

From equation (CIO) write the last derivative in equation (B4) as 

fel ■ frM < c27) 


Since riy = riy^TjP) , 


/dn v \ /dn v \ / ^ n v\ / St\ 

Ws = \5 p"" /T + \^)A^)s 


(C28) 


Now, (Sp/St) s has been evaluated in equations (C26) and (5). It can be 
expressed as (Sp/dT) s = ep/T and the other derivatives of equation (C27) 
are given by equations (C12) and (C2l) . With these substitutions equa- 
tion (C27) becomes 

/^n±\ fn v (AH v /RT) n v ] /Sp\ 

V5T/ S = [fell - \) ' P(1 - V J 


(C29) 
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Also, 


( &t\ _ /St\ /8p\ 

N = NslNs 

-£(£)l 


(C30) 


Substituting equations (C29) and (C30) into equation (B8) and solving 
for (dp/dp) s gives 




■{r 


Xv 


* * A;WJ * ‘I 


Eliminating e by means of equation (5) gives 

/Sp\ _ P 

Ws' 




r ! , x v j 

/ah v V 

'2 

\ 

o < 

1 

L r-Tx^ 

[rt ) 



H 

1 

\>f 

1 

i — 1 

, + \ 

(fj r 

) 2 ' 

’ > 



_ fr 

\RT- ( 

1 J 

j 


Finally 


_ p /& p\ 

'Hh 


i - 


1 + 


1 - X, 


Ml 


or 


6 fr + 


1 - X v 

The speed of sound in a region of phase change is then 


r 1 + ^ i 

1 

i - V 

)\ 


-Sr featY H 
1 - ^ V RT / J 


2 P 

c = gTp 


(C31) 


(C32) 


(C33) 


(C33a) 


(C34) 
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Equation (C33) was derived by assuming that the B2O3 present con- 
denses. If the mixture is superheated or if expansion is frozen, the 
terms X v /(l - X v ) and X^/R^l - X v ) in equation (C3l) both Danish, 
since they result only from considering condensation. Equation (9) in 
the DISCUSSION then replaces equation (C33) . 1 ' 


Biff s^ential Form of General Energy Equation 
For an adiabatic process the general energy equation is 

(U*V 2gj) + (h/Zn k m k ) = Constant 
This can be differentiated to give 

UdU 


gJ 


+ dh = 0 


For an isentropic process, equations (Bl) and (B4) give 

fbh\ RTZni 
VVs = ~ 

Combining equations (C36) and (C37) gives 

Zn 


(UdU) = &L 

s Zn k m k 



dP = -gJRT ^ 

s ^ n k m k ^ 


But if m m is defined as equation (C38) becomes 

(U®) s = - d (ln P ) 


(C35) 


(C36) 


(C37) 


(C38) 


(C39) 


Another form of this expression results from combining equation (C39) 
with equation (4) : v ' 


(UdU) s = - ^ dT 


(C40) 
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APPENDIX D 


CALCULATION OF THE EXAMPLE 

The calculation of the example consists of three parts. First, the 
combustor outlet conditions are determined. Second, the conditions in 
the nozzle at the theoretical saturation point are found. Third, the 
conditions through the remainder of the nozzle with equilibrium condensa- 
tion are determined by a series of stepwise calculations based upon the 
equations developed herein. 


Combustor Outlet Conditions 


The composition of the burned mixture, neglecting dissociation, is 
found as follows. The chemical formula of a fuel or mixture of fuels 
containing boron, carbon, and hydrogen can be represented as B a HpC\, 
where a, p, and X are integers or decimals. The atom fractions are 
then: for boron, a/ (a + p + X) ; for hydrogen, p/(a + p + X)j for carbon, 
X/(a + p + X) . These atomic ratios and the triangular chart of figure 2 
are used to find the moles ®2®3> ^gO ^ 2 ®’ and n C 02 ^ 9 2 

formed by burning the fuel stoichiometrically in 1 pound of air. The 
composition for all equivalence ratios less than stoichiometric is then 

moles BgOj/lb air 


$n B 


2°3 


$ 



moles H^O/lb air 


*n> 0 moles CO^/lb air 


0.00726(1-$) moles 02 /Lb air 
0.02740 moles Ng/lb air 


The stoichiometric fuel-air ratio used in computing $ can be found by 
means of the atom fraction and figure 3. 


For ethyldecaborane (B]_q c 2 h 18) > tile fuel considered in the example, 
the atom ratios are 0.333 for boron, 0.600 for hydrogen, and 0.067 for 
carbon. From figure 2 the stoichiometric mole numbers are found to be: 
n i 2 03 > °* 00259 > n H 2 0' °* 004:66 ; and n C0 2 ' 0, ° 0103 * ^' rom figure 3 the 


stoichiometric fuel-air ratio is found to be 0.0778. 
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If the combustion temperature is less than the saturation tempera- 
ture for B2O3 vapor, both condensed and gaseous B2O3 will be present. 

The criterion for the presence of the condensed phase is as follows: 

n B 2 03 P /HJlV (Dl) 


where 


Zn k = ^(n^ 2 o 3 + + n^) + 0.00726 (l - *) + 0.02740 (D2) 

The vapor pressure is plotted against temperature in figure 4. If in- 
equality (Dl) indicates no condensed material, the amount of B2O3 vapor 
is merely 


■V - *°b 2 o 3 


moles/lb air 


(D3) 


If the inequality indicates condensed material, the amount of B2O3 vapor 


n =Ev^L 

V P " P, 


m 


where 

Zn g = 2n k -4>r^ 2()3 

and the amount of condensed B2O3 is merely 

n c = * n B 2 0 3 ” “v 

The amount of noncondensing gases and B2O3 vapor is 

Za^ = Zn g + riy 


(D5) 


(D6) 


(D7) 


The total enthalpy of the burned mixture for a fuel containing boron, 
carbon, and hydrogen is then 

^ + = nc ^ H T^B 2 0 3 ,c + n v(H§)B 2 0 3 ,v + ^HgOt^HgO + + 


0.00726(1 -$)(H^)q + 0.02740(H^)jj Btu/lb air 

2 2 


(D8) 
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Total enthalpies of the constituents are plotted in figure 5. They in- 
clude chemical energy. Data for B 2 0 3 are from reference 6, while data 
for the other materials are from reference 1. 

Conservation of energy across the combustor is expressed as 



Enthalpy of air is given in figure 6. Enthalpies of some fuels of inter- 
est taken from reference 2 are as follows: 


Fuel 

Formula 

Phase 

Tempera- 

ture, 

°R 

Assigned 

enthalpy, 

*f, 

Btu/lb 

Boron 

B 

Crystal 


28,843 

Diborane 

b 2 h 6 

Gas 

536.7 

36,575 

Ethyldecaborane 

B 10 C 2 H 18 

Liquid 

536.7 

36,263 

Hydrocarbon 

fuel 

C x H 2x 

Liquid 


20,000 

Pentaborane 

b 5% 

Liquid 

536.7 

33,828 


Enthalpies of other fuels containing boron, carbon, and hydrogen in known 
proportions and having experimentally determined or empirically calculated 
heats of combustion can be found from the equation 



For heats of combustion determined at 298.16° K (536.7° R), (h£)q is 
24.642 Btu/lb mole, (H§) COs is 4.026 Btu/lb mole, (H§) B2 o 3 ( crystal) 
is 91.605 Btu/lb mole, and ( vapor ) is 24.642 Btu/lb mole. 

Combustion pressure P c and velocity U c required to evaluate ex- 
pressions (Dl) , (D3) , and (D9) are found by means of the momentum 
equation 



(DIO) 
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and the continuity equation 

^ ( X + P B% * Wc (^) c M 

where the static-pressure loss due to flameholder and combustor wall 
friction is represented by the ratio P-^/P-g. The mean molecular weight 
rr^ of the combustion products at station C can be expressed as 

m m ^ c = 69.64(X C + Xy) + 18.016 X^q + 44.01 Xqq^ + 32Xq^ + 28.016 X^ 

(D12) 


Also, (R/m a );g ^- s 0-06888 for the air composition assumed herein. 

The determination of the combustion temperature requires that equa- 
tions (D9) , (DIO) , and (Dll) be satisfied. Furthermore, if inequality 
(Dl) shows the presence of condensed B2O3, the proper amounts of n c 
and riy determined from equations (D4), (D5), and (D6) must be used in 
equation (D9) . The solution for a single combustion temperature may 
require several iterations. 

For the example being considered the combustor inlet temperature 
T-g was found to be 1590° R, while the inlet pressure was 4.747 atmos- 
pheres. By iteration of equations (D9), (DIO), and (Dll) the flame tem- 
perature Tq was found to be 4400° R neglecting dissociation. Velocity 
Uq was 980 feet per second, and pressure Pq was 4.524 atmospheres. 
Since inequality (Dl) revealed that no condensed B2O3 was present, the 
saturation condition for B2O3 may be reached in the expansion process. 


Conditions in Nozzle at Saturation 


The pressure and temperature at the saturation point can be found 
from conditions at the combustor outlet. For the combustion temperature 
T c and the parameter ^c/^ n k-’ tlle ra ‘ tio P]_/Pq is reacJ - from 

figure 7. The ratio T^/Tq can Be calculated from equation (13) using 


an e of 5.15. The subscript 1 refers to the condition where satura- 
tion is reached. For the example herein <£n^ Q ^c/^ n k '*' s 0-199. This 

value, together with the combustion temperature of 4400° R, determines a 
P]_/Pq of 0.637 and a T-^/Tq of 0.916. Thus P-^ and T^ are 2.885 
atmospheres and 4032° R, respectively. 


If the inequality (Dl) shows the presence of condensed B2O3 at the 
combustor outlet conditions, figure 7 cannot be used. The derivation 
of the equation from which figure 7 is constructed appears at the end of 
this appendix. 
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Expansion with Phase Equilibrium 

The equations developed in appendixes B and C are used to calculate 
conditions through the region of phase change in the nozzle. The detailed 
calculation procedure is illustrated in table II. The reader, by follow- 
ing the headings of table II, can make similar calculations for isentropic 
expansion with phase change. Alternate procedures are provided in table 
II for those cases where the mixture is superheated and those cases where 
the mixture is saturated. A test is provided for determining which case 
is present at any condition. 

At the precise point where the condensation can commence, the com- 
putation can be made either way. The temperature and pressure in rows 1 
and 1(a) of table II for the example are such that the test for condensed 
B 2 0 3 s k° ws "this to be the condition where saturation is reached. This is 
a consequence of having used figure 7 to determine the saturation point 
for the example. Row 1 was worked assuming that the test for condensed 
material showed the mixture to be very slightly superheated. Row l(a), 
for the same temperature and pressure as row 1, was worked with the 
assumption that the mixture was very slightly below the point of satura- 
tion. Column 15 of table II can be expressed as 

k( R P k = r[ X ^ (C p)b 2 0 3 ,v + X c^B 2 0 3;C + X H 2 0( C p)h 2 0 + 

x co 2 ( c p)co 2 + X 0 2 ( C p)o 2 + % 2 ( c p)n 2 ] ( D13 ) 

The variation of the constant-pressure specific heats Cp with tempera- 
ture for the products is given in figure 8. Inspection of column 15 in 
table II shows that summation (D13) could have been approximated by a 
value of about 5 over the entire range of conditions in the table. Values 
of AHy/RT for B 2 0 3 , column 16, are shown in figure 9 as a function of 
temperature . 


Calculation of Nozzle Contours 

Consider the case of constant change of velocity with Mach number. 
The distance from the nozzle inlet to any station rj for a nozzle 2 feet 
in length is 


\ ■ upfjJJ < u n - u c> 

For the case of constant change of Mach number with length there is 
a discontinuity in Mach number at the point of condensation. This is 
circumvented by calculating the distance l as follows 
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l 


■n 


S 

M D " M l(a) + M i ' Mg 


(Mtj - Me) 


Derivation of Equation Used to Construct Figure 7 

If inequality (Dl) shows no condensed B^O^ at the combustor outlet > 
the partial pressure of BgOj in the nozzle can be represented at any 
point down to the saturation condition by 

*°B203 P 

Zn k 

At the saturation point (station l) the partial pressure is equal to the 
partial pressure of the vapor in equilibrium with the condensed phase. 

The latter can be represented approximately over quite a range of temper- 
ature by 

e 14.87 - 68,300/T x 

Equating these expressions for the saturation point gives 

$n;B 2 0 3 Pl 14.87 - 68,300/T n 

n k 

Assume that composition does not change between station C and station 1. 
It is then possible to define a mean e^ r for this interval such that 

T c/ T i = ( p c/ p i) 1//Tfr 

Equation (D14) can then be rearranged as 



n k 



14.87 

e 



(D15) 


T 


fr 


A constant value for 


of 5.15 was chosen in constructing figure 7. 


28 


NACA RM E57C11 


REFERENCES 

1. Huff, Vearl N., Gordon, Sanford, and Morrell, Virginia E.: General 

Method and Thermodynamic Tables for Computation of Equilibrium Com- 
position and Temperature of Chemical Reactions. NACA Rep. 1037, 
1951. (Supersedes NACA TN's 2113 and 2161.) 

2. Hall, Eldon W., and Weber, Richard J. : Tables and Charts for Thermo- 

dynamic Calculations Involving Air and Fuels Containing Boron, 
Carbon, Hydrogen, and Oxygen. NACA RM E56B27, 1956. 

3. Maxwell, W. R., Dickinson, W. , and Caldin, E. F. : Adiabatic Expansion 

of a Gas Stream Containing Solid Particles. Aircraft Eng., vol. 
XVIII, no. 212, Oct. 1946, pp. 350-351. 

4. Glasstone, Samuel: Thermodynamics for Chemists. D. Van Nostrand Co., 

Inc., 1947. 

5. Setze, Paul C.: A Study of Liquid Boron Oxide Particle Growth Rates 

in a Gas Stream for a Simulated Jet Engine Combustor. NACA RM 
E55I20a, 1957. 

6. Setze, Paul C.: A Review of the Physical and Thermodynamic Properties 

of Boric Oxide. NACA RM E57B14, 1957. 


852 f 


4248 


MCA RM E57C11 


29 


TABLE I. - PROPERTIES OP EXPANDING COMBUSTION PRODUCTS AT POINTS 
WITHIN REGION OF BORIC OXIDE (B 2 0 3 ) CONDENSATION 


Nozzle 

Static 

Static 

Veloc- 

Speed 

Mach 

Coefficient 

Coefficient 

sta- 

pres- 

temper- 

ity. 

of 

number. 

of pressure- 

of pressure- 

tion. 

sure^ 

ature, 

U, 

sound. 

M 

density 

temperature 

n 

p, 

T, ' 

ft/sec 

C; 


variation. 

variation. 


atm 

°R 


ft/sec 


r 

6 

1 

2.885 

4032 

2705 

2910 

0.9296 

1.253 

4.947 

1(a) 

2.885 

4032 

2705 

2787 

.9706 

1.149 

10.048 

2 

1.817 

3851 

3634 

2709 

1.341 

1.151 

9.584 

3 

1.144 

3670 

4325 

2636 

1.641 

1.156 

8.978 

4 

.720 

3486 

4880 

2568 

1.900 

1.166 

8.090 

5 

.454 

3293 

5341 

2507 

2.130 

1.183 

7.045 

6 

.286 

3084 

5738 

2445 

2.347 

1.207 

6.071 

7 

.180 

2858 

6090 

2378 

2.565 

1.230 

5.408 

8 

.113 

2622 

6409 

2287 

2.802 

1.244 

5.115 

a 9 

.0712 

2396 

6689 

2191 

3.053 

1.250 

4.993 


a Nozzle outlet, station D. 


TABLE II. - CALCULATION PROCEDURE FOR ISENTROPIC EXPANSION PROCESS WITH EQUILIBRIUM CONDENSATION 


Nozzle 

sta- 

tion, 

7? 

1 

2 

3 

4 

5 1 

1 5 1 


8 

5_J 

1 2 

9 

10 

P, 

atm 


T, 

1 

tPn B 2 03 

Zn k 

Pv 

(fig. 4), 
atm 

All 

to 

H 

if (3) 

1 (4) 

x v , 

(®)/(8) 

>- X 

\o 

CD 

P " P v . 
(l)-(4) 

U)zn g 

n o- 

* n B 2 ° 3 

-(6) 

Zn^, 

(6)+In g 

•V 

* n B 2 0 3 

Znj, 

(6)+Zn g 

(2)77-1 

( 1 )77 

(19)77-1 

» 

^7?-l 

°R 

“W 

1 

2.885 


4032 


0.127 

0.127 





0.00155 

0.03527 

0.0439 

0 

1(a) 

2.885 


4032 


.127 

.127 

2.758 

0.00155 

0 

0.03527 



.0439 

0 

2 

1.817 


3851 


.080 

.058 

1.759 

.00111 

.00044 

.03483 



.0319 

.0126 

3 

1.144 


3670 


.050 

.0245 

1.1195 

.00074 

.00081 

.03446 



.0215 

.0235 

4 

.720 


3486 


.0316 

.00925 

.71075 

.00044 

.00111 

.03416 



.0129 

.0325 

5 

.454 


3293 


.0199 

.00289 

.45111 

.00022 

.00133 

.03394 



.0065 

.0392 

6 

.286 


3084 


.0126 

.000700 

.28530 

.000083 

.00147 

.03380 



.0025 

.0435 

7 

.180 


2858 


.0079 

.000115 

.17988 

.000022 

.00153 

.03374 



.000652 

.0453 

6 

.113 


2622 


.0050 

.0000124 

.11299 

.0000037 

.001546 

.033724 



.000110 

.0458 

9 

.0712 


2396 


.00312 

.00000097 

.071199 

.00000046 

.0015495 

.0337205 



.00001 

.0460 


Nozzle 

sta- 

tion, 

n 


1(a) 

2 

3 

4 

5 

6 

7 

8 
9 


Nozzle 

sta- 

tion, 

7? 


1 

1(a) 

2 

3 

4 

5 

6 

7 

8 
9 


x h 2 o> 


4n, 


H 2 0 

X sr 


l C0 2 ' 
* n C0 2 


0.0794 

.0794 

.0803 

.0811 

.0819 

.0824 

.0827 

.0829 

.0829 

.0829 


23 


2X k m k 


29.66 

29.66 

30.03 

30.36 

30.63 

30.82 

30.96 

31.01 

31.01 

31.02 


12 


W 


0.0176 

.0176 

.0177 

.0179 

.0181 

.0182 

.0183 

.0183 

.0183 

.0183 


13 


V 

1-4 

T*r 


0.0822 

.0822 

.0834 

.0843 

.0850 

.0856 

.0859 

.0861 

.0861 

.0861 


14 


x n 2 > 

0.02740 

T*J 


0.7768 

.7768 

.7867 

.7951 

.8021 

.8073 

.8106 

.8121 

.8125 

.8126 


15 


*X k Cg, k 

~R 


4.947 

4.947 

5.020 

5.077 

5.121 

5.139 

5.135 

5.101 

5.050 

4.985 


24 


lb/cu ft 


0.02906 

.02906 

.01940 

.01296 

.008664 

.005819 

.003932 

.002675 

.001830 

.001262 


25 


( 19 )t?-1 + ( 19 )t? 
5 


9.8160 

9.281 

8.534 

7.5675 

6.558 

5.7395 

5.2615 

5.054 


16 


1 iQ 


19 


If (3) > (4) 


AH V 

RT 

(fig- 9) 


16.7 

17.5 

18.55 

19.65 

20.95 
22.60 
24.60 
27.10 

29.95 


1 - x v , 

1 - (9) 


0.9561 

.9681 

.9785 

.9871 

.9935 

.9975 

.99935 

.99989 

.99999 


26 


AT, 

(2) 7 ?_i-(2) 7? 


181 

181 

184 

193 

209 

226 

236 

226 


(9)(16) 


T TV 


0.7668 

.5766 

.4076 

.2568 

.1371 

.0566 

.0160 

.00298 

.00030 


(15)+(16)(18) 
1 + (18) 


10.048 

9.584 

8.978 

8.090 

7.045 

6.071 

5.408 

5.115 

4.993 


27 


99 ,492(25 ) ( 26 ) 
(23) 


5.8879x10® 

5.5024 

5.1018 

4.7143 

4.4044 

4.1611 

3.9833 

3.6629 


If (3)<(4) 


(a) 


28 


(27)ff+(28)jr}_i 


7.31769X10® 

7.31769 

13.2056 

18.7080 

23.8098 

28.5241 

32.9285 

37.0896 

41.0729 

44.7358 


20 


22 


If (3) > (4) 


1 

1/(17) 


1.0459 

1.0330 

1.0220 

1.0130 

1.0065 

1.0025 

1.0006 

1.0001 

1.0 


29 


U, 

V(28), 

ft/sec 


2705.1 

2705.1 

3634.0 

4325.3 

4879.5 

5340.8 

5738 

6090 

6409 

6689 




0.1758 

.1645 

.1568 

.1554 

.1614 

.1740 

.1879 

.1961 

.2003 


30 


A, 

1 

(24) (29) " 
sq ft 
(lb/sec) 


0.01272 

.01272 

.01418 

.01784 

.02365 

.03218 

.04433 

.06139 

.08525 

.11846 


y> 

(207- C2l) 


1.149 

1.151 

1.156 

1.166 

1.183 

1.207 

1.230 

1.244 

1.250 


22 


If (3)<(4) 


TT^T 


1.253 


31 


v ra fgp r > 


(24) 

ft 

sec 


2910 

2787 

2709 

2636 

2568 

2507 

2445 

2374 

2287 

2191 


32 


M, 

(29) 

(31) 


0.9296 

.9706 

1.341 

1.641 

1.900 

2.130 

2.347 

2.565 

2.802 

3.053 


OJ 

o 


a Values same as in column 15. 
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Figure 1. - Nozzle contour for constant change of velocity with distance and for constant change 
of Mach number with distance. Ramjet engine at 60,000 feet and flight Mach number of 4.0; 

ethyldecaborane fuel; equivalence ratio, 0.6. w 
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Figure 2. - Moles of boric oxide, carbon dioxide, and water formed by stoichiometric burning of fuel 
in 1 pound of air. 


QfZf 


NACA EM E57C11 


4248 





.50 



Figure 3. - Stoichiometric fuel-air ratio for fuel containing boron, carbon, and hydrogen. 
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(a) Boric oxide. 
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Figure 5. - Variation of total enthalpy per mole with temperature for combustion-product 
constituents. * 
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(b) Water, carbon dioxide, oxygen, and nitrogen. 
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Figure 5. - Concluded. Variation of total enthalpy per mole with temperature for combustion- product constituents. 
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Temperature, °R 

Figure 6. - Variation of total enthalpy of air with temperature. 
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Figure 7 . - Chart for determining conditions at saturation temperature 
from known combustor outlet conditions. 
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(a) Boric oxide. 
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Figure 8. - Variation of molar constant-pressure specific heat with temperature for combustion- 
product constituents. 
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(b) Carbon dioxide, water, oxygen, and nitrogen. 


Figure 8. - Concluded. Variation of molar constant- pressure specific heat with temperature for combust ion- product 
constituents . 
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Figure 9. 


Variation of AH V /RT 


with temperature for 
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